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EXECUTIVE  SUMMARY 


A  microwave  landing  system  (MLS)  mathematical  modeling  study  was  performed  for 
Runway  22L,  Midway  Airport,  Chicago,  Illinois.  A  brief  model  description  is 
included  in  the  report  explaining  its  organization  and  capabilites.  This  study 
considered  the  static  effects  of  fixed  objects  and  the  transient  effects  of 
aircraft  likely  to  be  found  at  Midway. 

Ten  buildings  and  10  aircraft  were  modeled  for  multipath  effects.  Six  buildings 
and  one  aircraft  were  modeled  for  shadowing  effects.  Results  indicate  that  MLS 
performance  would  be  within  error  tolerances.  However,  plots  of  the  orbital 
flightpath  indicate  the  possibility  of  an  out-of-tolerance  condition  between  -24° 
and  -26°  due  to  the  shadowing  effects  of  buildings  7  and  8  ^new  Beckett  Aviation 
hangar).  These  effects  will  be  negligible  if  the  approach  procedure  keeps  the 
aircraft  at  3000  feet  mean  sea  level  (m.s.l.)  or  higher  in  the  vicinity  of  these 
radials. 


INTRODUCTION 


OBJECTIVE. 

To  identify  the  magnitude  of  the  potential  derogatory  effects  of  the  Midway 
Airport  environment  at  Chicago,  Illinois,  upon  a  microwave  landing  system  (MLS) 
precision  approach  to  runway  22L. 

BACKGROUND. 


Precision  MLS  approach  guidance  path  performance  may  be  derogated  by  the  effects 
of  reflections  (multipath)  from  buildings,  aircraft,  and  ground,  as  well  as  the 
diffraction  and  blockage  effects  (shadowing)  of  buildings  and  aircraft.  A 
computer  program  (mathematical  model)  has  been  developed  to  simulate  these 
effects  based  upon  user  inputs  describing  the  applicable  airport  scenario.  This 
program  is  described  by  the  reports  listed  in  the  bibliography.  Appendix  A 
provides  a  summary  description  of  the  math  model. 

Although  Congress  has  mandated  an  instrument  landing  system  (ILS)  installation 
for  Midway  Airport  runway  22L  in  the  FY-87  supplemental  budget,  the  Great  Lakes 
region  has  determined  that  an  ILS  installation  is  impractical  to  site.  This 
determination  was  based  upon  a  hostile  multipath  environment  and  a  requited  3.6° 
glidepath  for  obstacle  clearance  of  the  downtown  Sears  building.  Subsequently, 
since  Midway  22L  already  exists  as  part  of  the  Hazeltine  MLS  contract,  it  has 
been  moved  up  in  priority  to  provide  an  MLS  approach  by 

November  1988.  This  mathematical  modeling  study  was  performed  to  determine  if 
there  are  any  significant  multipath  effects  which  would  prevent  the  commissioning 
of  a  runway  22L  MLS  installation. 


MLS  MODEL  INPUT  CONSIDERATIONS 


This  simulation  was  based  upon  a  "quick-look"  philosophy  in  that  several 
short-cuts  were  applied  in  determining  input  data.  The  sites  selected  were  based 
on  nominal  locations  with  only  a  cursory  consideration  of  all  siting  regulations 
involved.  The  coordinates  used  for  buildings  and  aircraft  were  from  the 
digitizing  tablet  supported  by  nominal  field  survey  measurements  and  will  differ 
from  values  obtained  by  a  rigorous  site  survey.  Building  heights  were  calculated 
from  the  field  survev  measurements,  although  the  reference  locations  were  not 
surveyed  exactly.  However,  this  is  not  expected  to  have  any  noticeable  impact  on 
the  output. 

1.  Coordinate  Systems.  The  MLS  mathematical  model  uses  the  runwav  centerline  as 
the  X-axis  with  the  "0"  value  of  the  X  and  Y  axes  corresponding  to  the  stop  end 
of  the  runwav.  The  Z-axis  "1"  reference  is  thosen  as  the  lowest  mean  sea  level 
•m.s.l.)  value  along  the  runwav  which  is  the  threshold  (displaced1  of  runwav  32L 
at  Midway.  MLS  coverage  requirements  and  error  plots  are  based  on  a  coordinate 
system  centered  on  the  MLS  datum  point  (point  on  runway  centerline  opposite  the 
Elevation  (EL)  subsvstem). 


2.  Approach  Paths.  The  approach  path  descent  is  simulated  beginning  at 
6.242  nautical  miles  t nmi )  from  threshold  (3000  feet,  m.s.l.),  continuing  along 
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the  3.6°  glide  slope  to  a  point  8  feet  above  the  runway  surface  and  then  over  the 
runway  to  the  stop  end  remaining  at  8  feet  above  the  runway  surface.  In 
addition,  the  lower  coverage  limit  is  checked  with  a  10  nmi  orbit  at  0.9° 
elevation  angle  referenced  to  the  MLS  datum  point  from  -40°  to  +40°  horizontally 
referenced  to  the  MLS  datum  point  and  the  runway  centerline. 

3.  Multipath  Effects.  A  total  of  11  buildings  was  simulated  for  multipath. 

The  10  existing  and  proposed  buildings  with  the  highest  multipath  levels  were 
included  in  the  final  model  runs.  The  buildings  modeled  with  corresponding 
reference  numbers  (used  in  figure  1  and  tables  2,  5,  and  6)  are  as  follows: 


Description  Based  on  Future  Airport  Layout  Plan 

Monarch  hangar 
Monarch  hangar 
Proposed  corporate  hangar 
Esmark  hangar  installation 
Monarch  and  Butler  hangars 
Midway  Airlines  hangar 
Beckett  new  hangar  office  area 
Beckett  new  hangar 
ATC  tower  north  side 
ATC  tower  west  side 


A  total  of  15  aircraft  (B-727's)  were  also  simulated  at  various  locations  on  the 
airport  taxiways  and  ramp  areas.  This  number  was  reduced  to  the  10  most  likely 
to  cause  problems  based  upon  multipath  levels.  Aircraft  locations  simulated  on 
the  ramp  areas  were  based  upon  aerial  photographs.  The  modeled  aircraft  with 
corresponding  reference  numbers  (used  in  figure  1  and  tables  2,  5,  and  6)  are  as 
follows : 


Description  Based  on  Future  Airport  Layout  Plan 

On  north  ramp  between  runways  22R  and  22L 
Holding  on  taxiwav  east  of  runwav  4R  stop  end 
Holding  on  taxiway  near  existing  displaced  runway  22L 
threshold 

Parked  north  of  concourse  C 
Parked  south  of  concourse  C 
Parked  north  of  concourse  B  near  tower 
On  runway  4R/22L  taxiway  just  south  of  runway  13R/31L 
taxiwav 

On  runwav  13R  just  north  of  runwav  4R/22L 
Parked  at  northwest  corner  of  concourse  B 
On  4R/22L  taxiway  just  north  of  south  taxiwav 


Although  the  downtown  Sears  building  and  other  nearby  skyscrapers  are  an 
obstruction  problem,  they  were  not  considered  in  this  modeling  studv  since  the 
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approach  path  modeled  did  not  extend  beyond  downtown  Chicago.  The  approach  path 
modeled  was  the  only  specific  path  available  since  an  approach  procedure  fox  a 
Midway  22L  MLS  approach  has  not  been  developed  yet.  When  the  approach  procedure 
is  determined,  additional  modeling  should  be  performed  for  that  approach  path. 

4.  Shadowing  Effects.  Shadowing  effects  include  both  blockage  by  an  object  and 
diffraction  around  the  object.  The  effects  of  a  shadowing  aircraft  on  the 
azimuth  subsystem  were  simulated  by  an  aircraft  taking  off  from  runway  13R.  To 
simulate  a  worst  case  condition,  the  interfering  aircraft  just  crossed  runway 
22L  when  the  aircraft  making  an  MLS  approach  to  runway  22L  was  over  the 
threshold.  The  effects  of  six  shadowing  buildings  were  also  simulated.  The 
buildings  included  were  identified  above  as  numbers  5,  6,  7.  8.  9.  and  10. 


DATA  PRESENTATION  AND  ANALYSIS 


Input  data  unique  to  this  scenario  are  listed  in  table  1,  MLS  Input  Parameters, 
and  table  2,  Scenario  Obstacle  Data.  These  tables  list  the  transmitter 
locations,  building  locations,  etc.,  used  for  this  simulation.  A  detailed 
explanation  of  the  various  input  parameters  is  provided  in  appendix  B.  An 
airport  scenario  map  showing  the  relationship  of  the  runway  and  transmitters  to 
the  multipath  sources  and  shadowing  objects  is  shown  in  figure  1.  Building 
locations  are  represented  by  rectangles  (wide  lines)  and  are  referenced  to 
table  2  by  the  adjacent  numbers.  Aircraft  locations  are  indicated  by  arrow 
shaped  symbols  and  may  be  referenced  to  table  2  by  the  number  following  the  "A." 
The  tip  of  the  arrow  indicates  the  aircraft  nose.  Waypoint  and  segment 
parameters  used  for  the  centerline  approach  are  listed  in  table  3.  Current lv. 
the  orbital  flightpath  can  only  be  simulated  by  a  series  of  segments.  The 
waypoint  and  segment  parameters  for  the  orbital  flightpath  are  provided  in 
table  4.  The  orbital  flightpath  was  simulated  at  10  nmi.  as  opposed  to  the 
20  nmi  coverage  limit,  in  order  to  minimize  computer  run  time  and  the  number  of 
data  points.  The  altitude  used  for  the  orbital  flightpath  places  the  aircraft 
at  the  lower  coverage  limit. 

The  maximum  values  for  multipath  from  the  ground,  buildings,  and  aircraft  are 
ranked  and  listed  in  table  3  for  the  centerline  approach  and  m  table  b  for  the 
orbital  flightpath.  Diagnostic  plots  show  the  "Multipath/Direct"  iM/D1  ratio  for 
all  Azimuth  (AZ),  EL,  Precision  Distance  Measuring  Equipment  tDME/P)) 
subsvstems.  In  addition,  ''Separation  Angle”  plots  are  provided  for  the  angle 
equipment,  and  "Relative  Time  Delay"  plots  are  provided  for  the  DME,P.  In  the 
upper  right  corner  of  the  diagnostic  plots  is  a  legend  indicating  which 
multipath  sources  correspond  with  the  plot  symbols.  The  legend  list  is  for  the 
highest  six  multipath  sources  and  is  ranked  accordingly.  The  solid  line  on  the 
M/D  plots  is  used  to  connect  the  data  points  for  the  highest  ranked  multipath 
source.  Values  ire  plotted  on  the  separation  ingle  plots  onl"  when  the 
corresponding  Multipath/Direct  'M/D'  ratio  is  ibove  --*0  decibels  ‘dB1.  The  solid 
line  on  the  separation  angle  plots  is  used  to  connect  the  svmbols  where  tne 
multipath  exists  continuously  over  several  samples. 

Figures  2  and  3  are  the  M/D  and  separation  angle  plots,  respectively,  for  the  AZ 
subsystem.  These  plots  have  an  expanded  X-axis  to  show  the  multipath  effects 
near  threshold  more  clearly  since  no  multipath  effects  occur  on  the  approach 
beyond  1  nmi  from  threshold.  The  six  highest  multipath  sources  t tanked  1  through 
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6  in  table  5  for  A2)  are  shown  on  these  plots  as  identified  by  the  legend  in  the 
upper  right  corner.  Although  building  8  shows  a  multipath  level  within  1  dB  of 
the  direct  signal  in  the  touchdown  zone,  the  separation  angle  exceeds  10°  in  this 
area.  The  multipath  from  buildings  7  and  8  occurs  further  down  the  runway,  and 
again  the  multipath  is  associated  with  large  separation  angles.  Similarly,  for 
the  remaining  buildings  and  aircraft,  the  separation  angle  for  the  multipath 
sources  is  large  enough  (greater  than  2  beam  widths)  and  the  M/D  ratio  is  low 
enough  to  prevent  any  significant  problems.  Therefore,  no  significant  AZ  errors 
are  expected  due  to  multipath  from  aircraft  or  buildings. 

The  M/D  and  time  delay  diagnostic  plots  for  the  DME/P  subsystem  are  provided  in 
figures  4  and  5,  respectively.  These  plots  also  have  an  expanded  X-axis  to  show 
the  multipath  effects  near  threshold  more  clearly  since  no  multipath  effects 
occur  on  the  approach  beyond  1  nmi  from  threshold.  The  six  highest  multipath 
sou  .  ?s  (ranked  1  through  6  in  table  5  for  DME/P)  are  shown  on  these  plots,  as 
iden^_fied  by  the  legend  in  the  upper  right  corner.  Building  6  is  shown  to  have 
a  multipath  ievei  approaching  that  of  the  direct  signal.  However,  the  associate i 
time  delay  is  in  excess  of  700  nanoseconds  (ns)  which  should  eliminate  any 
errors.  The  multipath  from  building  8  is  shown  to  come  within  2  dB  of  the  direct 
signal  in  an  area  where  the  time  delay  is  about  300  ns.  This  could  cause  some 
accuracy  errors.  However,  this  would  occur  approximately  3000  feet  past 
threshold  and,  therefore,  should  not  be  a  concern.  For  the  remaining  buildings 
and  aircraft  no  accuracy  errors  are  expected  due  to  low  multipath  levels  'more 
than  3  dB  below  the  direct)  or  long  time  delays  (greater  than  350  ns). 

Diagnostic  plots  for  the  EL  subsystem  are  shown  by  figures  6  and  7.  These  plots 
also  have  an  expanded  X-axis  to  show  the  multipath  effects  near  threshold  more 
clearly  since  no  multipath  effects  occur  on  the  approach  beyond  1  nmi  from 
threshold.  The  six  highest  multipath  sources  t ranked  1  through  b  in  table  5  ioi 
EL/  are  snown  on  these  piots,  as  laentitied  by  the  legend  m  the  upper  ri0hc 
corner.  These  plots  must  be  examined  with  caution  to  avoid  anv  misconceptions. 
High  multipath  levels  are  seen  behind  the  EL  antenna  and  should  be  ignored 
because  the  model  is  omnidirectional  at  this  stage.  Antenna  directivity  is 
considered  subsequently  in  the  system  part  of  the  model.  The  highest  level  of 
multipath  to  consider  is  the  in-beam  multipath  from  aircraft  b.  Howevei .  the  low 
M/D  ratio  (below  -19  dB)  should  preclude  any  errors. 

Shadowing  effects  to  the  AZ  subsystem  are  shown  in  figure  8.  The  amplitude 
fluctuations  are  caused  by  ground  lobing  effects  since  the  shadowing  simulation 
includes  ground  reflection  computations.  The  above  comments  also  ipplv  ro  the 
DME/P  shadowing  plot  in  figure  ■) .  ,<u  snaaowing  effects  are  snown  on  the  ZL  pint 

in  figure  10  since  the  aircraft  are  behind  the  EL  subsystem  and  the  beam  is 
shaped  to  minimize  ground  reflections. 

The  Time  Reference  Scanning  Beam  vTRSB)  receiver  simulation  routine  processes  the 
multipath  generated  in  the  oroDagation  part  of  the  model  to  determine  AZ  and  EL 
angle  errors.  This  processing  is  modeled  after  an  actual  receiver  flow  c.iart  mu 
includes  such  features  as  dwell  gate  or  split  gate  processing  'uweil  gate  used 
for  this  simulation/,  acquisition,  tracking,  system,  flags,  coast  mode,  md  slew 
rate  limiting.  The  transmitter  radiation  patterns  Haze  It  me  2°  AZ  and  1 . V  FI 
and  aircraft  antenna  patterns  are  also  applied  at  this  point  m  the  processing. 
The  raw  error  from  the  simulation  is  further  processed  bv  passing  it  through  Path 
Following  Error  (PFE)  and  Control  Motion  Noise  (CMN)  filters.  The  equations  used 
to  implement  the  filters  are  based  on  the  application  of  a  Bilinear 
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Transformation  to  the  transfer  function  (see  appendix  C).  Other  than  the 
addition  of  the  10  radians  per  second  low-pass  filter  to  the  CMN  calculation, 
the  equations  are  equivalent  to  those  being  used  for  other  MLS  data  processing 
activities.  Simulation  of  a  DME/P  interrogator  is  not  currently  included  in  the 
model  program. 

Figures  11,  12,  and  13  are  the  Raw  error,  PFE  filtered,  and  CMN  filtered  plots, 
respectively,  for  the  AZ  subsystem.  Although  considerable  raw  errors  are 
observed  in  figure  11,  the  filtered  errors  are  well  within  the  tolerance  limits 
as  shown  in  figures  12  and  13.  The  EL  subsystem  Raw  error,  PFE,  and  CMN  plots 
are  provided  in  figures  14,  15,  and  16.  The  spike  near  the  end  of  the  valid 
elevation  information  is  attributed  to  conical  angle  effects.  Due  to  the  lack 
of  significant  multipath  effects,  no  out-of-tolerance  errors  are  evident. 

However,  the  errors  shown  are  attributed  to  the  air  traffic  control  ( ATC  >  tower 
and  building  4  (scheduled  for  demolition). 

To  check  for  potential  problems  in  the  coverage  area  away  from  centerline,  an 
orbit  was  simulated  at  the  lower  coverage  limit  (0.9°  964  feet  above  ground'. 

This  type  of  simulation  also  provides  the  diagnostic,  raw,  and  filtered  error 
plots.  Figures  17  and  18  show  the  orbital  M/D  and  separation  angle  plots  for  the 
AZ  subsystem.  The  azimuth  radials  which  would  be  subjected  to  multipath  can  be 
readily  identified  in  figure  17  as  -42°,  -27°  through  -20°,  and  *3°.  The 
highest  multipath  observed  occurs  at  -24.8°  from  building  8  (see  table  6  for  AZ  i 
and  exceeds  the  level  of  the  direct  signal  (3.23  dB>.  Building  7,  which  is  part 
of  building  8,  also  provides  a  multipath  level  in  excess  of  the  direct  signal 
(1.73  dB).  The  separation  angles  for  these  buildinrs  are  much  greater  than  2 
beam  widths,  and  no  out-of  tolerance  errors  are  expected.  The  multipath  levels 
from  the  other  buildings  are  low  enough  and  the  separation  angles  large  enough  to 
preclude  any  significant  errors.  Orbital  M/D  and  time  delay  plots  for  the  DME  P 
are  provided  by  figures  19  and  20.  Radials  where  multipath  could  affect  the 
DME/P  are  easily  seen  in  figure  19  to  be  the  same  as  the  AZ.  These  multipath 
levels  are  all  more  than  3  dB  below  the  direct  signal  which  should  preclude  non¬ 
significant  accuracy  effects.  However,  due  to  the  short  time  delays  near  -41“ 
from  buildings  1,  2,  and  3  (less  than  350  ns),  this  issue  should  be  addressed  bv 
the  contractor.  Building  4  is  scheduled  for  demolition;  therefore,  the  effects 
near  +4°  may  be  disregarded.  EL  subsystem  plots  for  M/D  and  separation  angle 
from  the  orbital  flightpath  are  shown  in  figures  21  and  22.  Potential  problems 
exist  from  the  ATC  tower  (building  9)  due  to  the  high  level  of  in-beam  multipath 
near  *12.7°.  However,  since  the  tower  is  located  about  82°  from  the  EL 
boresight,  the  amount  of  signal  available  at  this  angle  due  to  the  application  of 
the  antenna  patterns  is  minimal,  and  any  resulting  errors  are  expected  to  be 
within  tolerances. 

Shadowing  plots  for  the  orbital  flightpath  begin  with  figure  23  for  the  AZ 
subsystem.  The  constant  5  dB  bias  in  this  plot  is  attributed  to  the  ground 
lobing  associated  with  the  shadowing  simulations.  The  amplitude  change  between 
-31°  to  -25°  is  attributed  to  the  diffraction  effects  of  building  5.  The 
shadowing  effects  of  building  b  are  seen  Detween  -15°  to  -12°.  The  ATC  tower 
effects  on  the  AZ  signal  can  be  observed  at  about  *12°.  The  effects  of  buildings 
7  and  8  are  between  22°  through  32°.  DME/P  shadowing  effects  are  presented  bv 
figure  24.  Due  to  the  difference  m  frequencies,  the  bias  and  magnitude  of 
shadowing  effects  are  different.  However,  the  buildings  causing  AZ  shadowing 
effects  also  affect  the  DME/P  at  the  same  angles  due  to  collocation  of  the  sites. 
The  only  shadowing  effect  apparent  to  the  EL  subsvstem  is  the  constant  bias  due 
to  ground  lobing  (see  figure  25). 


Raw  error,  PFE  filtered,  and  CMN  filtered  plots  are  provided  for  the  AZ 
subsystem  in  figures  26,  27,  and  28.  The  errors  generated  near  -25°  are 
attributed  to  the  multipath  effects  of  building  8,  whereas,  the  effects  near  *25° 
are  caused  by  the  diffraction  effects  (shadowing)  of  building  8.  The 
shadowingeffects  of  the  ATC  tower  are  observed  near  +12°.  The  shadowing  effects 
of  building  8  cause  an  out-of-tolerance  condition  between  22°  and  30°  as  shown  bv 
figures  27  and  28.  Otherwise,  all  errors  are  within  tolerances.  Orbits 
simulated  at  higher  altitudes  of  2000  and  3000  feet  above  ground  showed 
considerably  less  shadowing  effects  and  resulted  in  AZ  errors  which  were  within 
tolerance  (the  plots  are  not  included).  The  Raw  error.  PFE,  and  CMN  plots 
generated  by  the  orbital  flightpath  for  the  EL  subsystem  are  shown  in  figures 
29,  30,  and  31,  respectively.  Errors  from  building  4  (which  is  scheduled  for 
demolition)  near  centerline  are  evident  on  the  EL  subsystem  plots.  However,  ill 
errors  fall  well  within  the  dashed  error  tolerance  lines  of  figures  30  and  31. 


SUMMARY 


The  approach  in  modeling  this  airport  was  to  identify  potential  probler  net.- 
prior  to  the  MLS  equipment  installation.  Subsystem  locations  modeled  aie  based 
upon  nominal  siting  considerations  and  are  expected  to  conform  with  the 
siting  criteria.  Building  coordinates  (X  and  Y>  were  obtained  fiom.  r  di^.’i-: 
tablet  with  additional  data  obtained  from  a  cursorv  site  survev.  Due  r  a; 
scale,  parallax,  and  tablet  resolution,  these  values  will  diffei  fir-  t r.  - 
obtained  by  a  rigorous  field  survey.  Building  heights  were  either  esr  t 
obstruction  chart  data  or  computed  from  field  survev  data.  Due  t,,  i-vimy  . 
made  in  the  model,  the  results  presented  here  are  considered  "uois^-  i- 
whereas,  the  actual  effects  of  multipath  and  shadowing  are  expe.  r-d  - 
than  shown.  However,  any  areas  of  the  approach  and  the  orbit  ;nd  r 
potential  problems  should  be  addressed  by  the  contractor  in  the  - :  • » 
report. 

In  this  scenario,  the  buildings  (modeled  as  smooth,  perfect  r**t>  • 
aircraft  create  mainly  out-of-beam  and  low  levels  of  multipath.  >hi. • 
location  is  expected  to  be  free  of  anv  significant  multipath  pit 
approach.  The  results  of  the  orbit  indicate  that  thp  AZ  and  DMF  c  ■-  r 
performance  could  be  affected  by  diffraction  « shadowing  >  from  builiing- 
(new  Beckett  Aviation  hangar).  An  out-of-tolerance  condition  i-'siii'-  . 
lower  coverage  limit  between  -24°  through  -2b°.  If  the  approach  p :  >  • 
keeps  the  aircraft  at  3000  feet  above  ground  or  higher  in  rhe  •  ic: t,  T  •  ■ 

radials,  however,  the  shadowing  effects  of  these  buildings  will  be  i-x.. 

We  recommend  further  modeling  of  the  approach  procedures  when  the\  ue  :  i , .  t . 


.••■.WV.-.S-.N-.V-.'.V. 
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TABLE  2.  MIDWAY  RUNWAY  22L,  SCENARIO  OBSTACLE  DATA 


TABLE  3.  MIDWAY  RUNWAY  22L,  CENTERLINE  APPROACH  FLIGHTPATH 
WAYPOINTS  AND  DATA 


FLIGHTPATH  TYPE:  SEGMENTED 
DATUM  COORDINATES: 

X:  46  IS.  Y:  0.  2: 


TABLE  OF  FLIGHTPATH  AND  WAYPOINT  DATA 


AYPT 

X -COORD 

Y-COORD 

2— COORD 

VELOCITY 

SAMPLING 

DISTANCE 

(NM) 

ID 

<FT) 

(FT) 

(FT) 

(FT/SEC) 

I NCR  (FT) 

ALONG  FP 

FROM  TH 

1 

43270.51 

0.00 

2441.01 

200 . 00 

40 . 00 

7.  13 

-0.  88 

2 

4626.  1 1 

0.  00 

9.71 

200 . 00 

40 . 00 

13.51 

-7.  25 

3 

3750.00 

0.00 

11.00 

200. 00 

40.00 

13.65 

-7.40 

4 

2950.00 

0.00 

13.00 

200. 00 

40.00 

13.78 

—7 . 53 

5 

1375.00 

0.00 

18.00 

200.00 

40.00 

14.04 

-7.79 

6 

550.00 

0.00 

22.00 

200.00 

40 . 00 

14.  19 

-7.92 

7 

0.00 

0.00 

22.00 

200.00 

40 . 00 

14.27 

-8.01 
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TABLE  6.  Ill  1)U AY  RUNWAY  221.,  MULT  1  PATH  RANKINGS  FOR  ORBITAL  FLICHTPATH 
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APPENDIX  A 


MICROWAVE  LANDING  SYSTEM  (MLS)  MATHEMATICAL  MODEL  DESCRIPTION 


The  MLS  mathematical  model  simulation  program  is  written  in  the  FORTRAN  77 
(ANSI  X3. 9-1978)  computer  language  and  has  successfully  been  used  on  computers  in 
the  United  States,  United  Kingdom,  and  the  Federal  Republic  of  Germany.  An  MLS 
simulation  may  be  considered  as  consisting  of  three  processes. 

1.  The  first  process  is  the  creation  of  a  formatted  input  file  which  defines  the 
airport  environment.  This  input  data  specifies  the  locations  and  composition  of 
reflecting  and  shadowing  obstacles,  terrain  features,  antenna  locations,  and  the 
simulated  flightpath. 

2.  The  second  process  is  known  as  the  propagation  model.  This  program 
determines  the  signals  at  the  receiver  for  each  point  along  the  flightpath, 
taking  into  account  the  various  multipath  reflections.  The  diagnostic  plots 
obtained  from  the  propagation  model  show  which  obstacles  could  cause  significant 
multipath  effects. 

3.  The  third  process  is  the  Time  Reference  Scanning  Beam  (TRSB)  system  and 
receiver  model.  This  part  of  the  simulation  computes  the  receiver  error  caused 
by  multipath  for  the  specified  ground  equipment  antenna  patterns,  aircraft 
antenna  pattern,  scan  format,  and  receiver  processing  algorithm.  The  raw  errors 
computed  by  the  system  model  are  passed  through  PFE  and  CMN  filter  algorithms  and 
plotted  along  with  the  applicable  tolerance  limits. 

Figure  A-l  shows  the  interrelationship  of  these  processes  to  the  total  MLS 
simulation.  (See  Bibliography  in  text  for  detailed  theory  and  description.! 

AIRPORT/ FLIGHTPATH  MODEL. 

The  formatted  input  file  consists  of  data  specified  for  the  particular  airport 
environment  being  considered.  At  the  Federal  Aviation  Administration  (FAA) 
Technical  Center,  this  information  is  currently  entered  in  an  interactive  session 
which  edits  a  standard  input  file  template  to  add  the  appropriate  data  for  input 
to  the  propagation  model.  The  degree  of  approximation  to  an  actual  airport 
environment  will  depend  heavily  upon  the  simulation  and  the  scatterer  geometrv. 
For  example,  in  a  case  where  the  multipath  is  out  of  beam  and  of  short  duration, 
hangars  might  be  represented  by  a  single  plate;  however,  to  closely  predict 
actual  system  performance  in  a  critical  multipath  situation,  it  would  be 
necessary  to  input  the  same  hangar  with  many  plates  representing  the  various 
electrical  properties  of  the  different  parts  of  the  hangar. 

PROPAGATION  MODEL. 
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Propagation  modeling  consists  of  executing  the  propagation  program  with  the 
formatted  input  file  as  input  data.  This  model  determines  the  multipath 
characteristics  of  the  specified  airport  environment.  The  numerical  results  from 
the  propagation  model  define  the  direct  signal;  signals  reflected  from  or 
diffracted  by  terrain,  buildings,  and  aircraft;  and  the  changes  in  the  direct 
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signal  characteristics  lue  to  shadowing  bv  runwav  humps.  buildings,  and  aircraft. 

The  type  of  multipath  considered  in  a  simulation  is  dependent  upon  and  determined 
by  the  input  parameters  specified.  Numerical  results  are  saved  in  two  different 
files:  one  file  for  further  processing  by  the  system  model  (data  set  S),  and  the 
other  for  plotting  of  the  multipath  data  (data  set  16).  The  propagation  model 
can  accommodate  the  multipath  types  specified  below. 

1.  Terrain  reflection  modeling.  The  terrain  is  typically  represented  by  a 
collection  of  rectangular  and  triangular  plates,  each  with  prescribed 
orientation,  roughness,  and  dielectric  constant.  By  varying  these  parameters, 
one  can  assess  the  sensitivity  of  performance  to  terrain  type  (e.g.,  dry  ground 
versus  snow).  The  multipath  levels  are  computed  either  by  a  numerical 
Kirchof f-Fresnel  integral  or  a  simplified  approximation. 

2.  Building  reflection  modeling.  Buildings  are  represented  by  one  or  more 
rectangular  plates  of  prescribed  orientation  and  surface  material.  The  various 
plates  represent  salient  features  of  a  building  such  as  the  doors  of  a  hangar. 

By  allowing  each  plate  to  have  a  different  surface  material  characterization, 
inhomogeneous  surfaces  (e.g.,  concrete  walls  with  metal  doors)  can  be  modeled. 
Consideration  is  also  made  for  secondary  ground  reflection  paths.  The  levels  are 
computed  assuming  Fresnel  diffraction  and  using  closed  form  Fresnel  integral 
expressions . 

3.  Aircraft  reflection  modeling.  For  aircraft,  it  is  essential  to  consider  the 
curvature  of  the  surfaces  as  this  tends  to  spread  the  reflections  over  a  much 
greater  region  than  would  be  the  case  with  flat  plates.  The  fuselages  and  tail 
fins  are  both  modeled  as  cylinders  or  a  section  thereof.  The  resulting  multipath 
levels  are  computed  bv  a  combination  of  Fresnel  diffraction  (integrals1  and 
geometric  optics. 

4.  Shadowing.  Shadowing  by  buildings  or  aircraft  causes  both  an  attenuation  and 
distortion  of  the  transmitted  wavefront.  Both  of  these  factors  are  consideied  m 
the  models  for  shadowing.  The  shadowing  obstacles  are  represented  by  one  or  more 
rectangular  plates  which  approximate  the  object  silhouette.  Similar  techniques 
have  been  successfully  used  in  studying  the  effects  of  widebodv  aircraft  on  the 
Instrument  Landing  System  iILS).  The  shadowing  of  the  azimuth  signal  bv  runway 
humps  requires  explicit  consideration  of  the  surface  curvature  and  is  computed  bv 
mathematical  algorithms  similar  to  those  of  aircraft  reflection  modeling. 

PROPAGATION  MODEL  DIAGNOSTIC  PLOTS.  Graphical  interface  subroutines  have  been 
developed  to  display  the  model  output  on  several  alternate  graphical  display 
devices.  The  graphical  routines  generally  used  to  display  the  propagation  model 
results  at  the  FAA  Technical  Center  are  from  the  TEKTRONIX  Plot  10  Terminal 
Control  System  and  CALCOMP  Preview  software.  These  routines  provide  easy  access 
to  the  gnohic  capabilities  of  the  Tektronix  4010  tvpe  Direct  View  Storage  Tube 
iDVST)  terminals.  Information  displayed  on  the  storage  tube  :ra\  be  copied  as 
desired  (via  a  hard  copy  unit)  to  provide  a  permanent  record  of  the  resuits. 
Graphical  output  from  the  multipath  model  consists  of  a  listing  of  the  input 
parameters  used  in  the  simulation,  the  flightpath  of  the  receiver,  an  .input’.  :ap 
showing  the  location  of  the  transmitters  and  obstacles,  and  multipath 
diagnostics.  These  diagnostics  display  relative  azimuth  (AZ),  distance 
measuring  equipment  (DME),  and  elevation  (EL)  multipath, direct  iM/D*  amplitude 


ratios  (for  the  maximum  component  of  the  several  multipath  components  from  a 
given  obstacle)  and  separation  angles  (time  delay  for  DME)  along  the  flightpath 
for  the  obstacles  generating  significant  multipath  components,  and  the  variation 
in  the  direct  signal  A Z,  DME,  or  EL  level  where  shadowing  is  involved. 

SYSTEM  MODEL. 

The  TRSB  system  and  the  MLS  receiver  algorithms  are  simulated  in  the  system 
model.  This  model  considers  the  received  signal  as  a  superposition  of  the 
received  direct  path  signal  and  a  number  of  replicas  (multipath)  of  it,  each 
having  its  own  amplitude,  delay,  angle,  and  Doppler  shift.  The  system  model  then 
determines  the  receiver  error  by  taking  into  account  the  nature  of  the 
transmitted  signals  and  the  antenna  patterns.  The  functional  form  of  the  beam 
waveform  is  determined  from  measured  or  theoretical  patterns  and  is  included  in 
the  model  as  a  function  subprogram.  By  superimposing  the  beam  patterns 
corresponding  to  the  various  signal  paths,  the  net  received  envelope  is 
determined. 

The  remainder  of  the  system  model  parallels  the  processing  by  the  receiver 
microprocessor.  A  tracking  gate  (dwell  gate  or  split  gate  can  be  simulated)  is 
centered  on  the  largest  consistent  envelope  peak  with  the  beam  arrival  angle 
derived  by  finding  the  times  at  which  the  leading  and  trailing  edges  of  the 
received  envelope  cross  a  threshold.  Various  checks  and  tracking  algorithms  are 
applied  to  each  measurement  before  it  is  presented  as  angle  data.  DME  is  not  a 
part  of  the  system  model  at  this  time. 

SYSTEM/RECEIVER  MODEL  PLOTTING. 

The  output  of  the  system  model  is  generally  displayed  on  a  Tektronix  DVST  using 
the  TEKTRONIX  Plot  10  and  CALCOMP  Preview  graphics  subroutines.  A  specific 
transmitter  (AZ  or  EL)  is  selected  for  plotting  and  the  raw  vand  static  If 
desired)  errors  generated  by  the  model  are  plotted  versus  the  distance  along  the 
flightpath.  The  raw  errors  may  be  processed  with  digital  PFE  and  CMN  filtering 
algorithms.  The  resultant  errors  may  be  plotted  against  the  applicable 
tolerance  limits.  An  analysis  of  the  filtered  error  plots  shows  whether  the 
system  is  suitable  for  commissioning. 


APPENDIX  B 


DESCRIPTION  OF  INPUT  PARAMETERS  DISPLAYED  IN  TABLES  1  AND  2 


PARAMETERS  REQUIRED  FOR  MULTIPATH  COMPUTATIONS. 

The  following  parameters  are  required  to  specify  the  airport  model  which  is 
employed  in  the  multipath  computation  section  of  the  program.  A  standard 
rectangular  coordinate  system  is  used,  where  the  X,Y-plane  is  in  the  plane  of 
the  runway,  the  X-axis  is  coincident  with  the  runway  center  line  and  the  Z-axis 
passes  through  the  stop  end  of  the  runway.  All  lengths,  frequencies,  and  times 
are  given  in  feet,  hertz  (Hz),  and  seconds,  respectively. 

♦TRANSMITTER  PARAMETERS  (AZIMUTH,  PRECISION  DISTANCE  MEASURING  EQUIPMENT,  AND 

elevatim). - 

1.  Azimuth  -  (X,Y,Z),  Elevation  -  (X.Y.Z),  DME/P  -  (X,Y,Z):  X,Y,Z  - 
coordinates  of  location  of  transmitter. 

2.  (AZ,  EL,  DME/P)  FREQUENCY:  Frequency  of  transmitter. 

♦PARAMETERS  USED  IN  COMPUTATION  OF  MULTIPATH  REFLECTIONS. 

♦Rectangular  Surface  Elements  are: 

1.  ID:  Two-character  surface  identification. 

2.  (XI,  Y1 ,  Z 1 ) ,  (X2,  Y2,  Z2),  (X3,  Y3,  Z3):  X,Y,Z-coordinates  of  two  corners, 
plus  X,Z  coordinates  of  third  corner,  in  increasing  order  of  magnitude  for  the 
X-coordinate,  for  each  rectangular  surface  element. 

3.  ERSR,  ERSI,  SH2S:  The  real  and  imaginary  relative  dielectric  constants,  and 
the  root-mean-square  roughness  height,  respecti vely ,  for  each  rectangular  surface 
element. 

4.  NRSPEC:  Rectangular  ground  surface  flag. 

♦Triangular  Surface  Elements  are: 

1.  ID:  Two-character  surface  identification. 

2.  (XI,  Y1,  Z1),  (X2,  Y2,  Z2),  (X3,  Y3,  13):  X,Y,Z-coordinates  of  the  three 
corners  of  each  triangular  surface  element,  in  increasing  order  of  magnitude  of 
the  X-coordinate. 


3.  ERSR,  ERSI,  SH2S:  The  real  and  imaginary  relative  dielectric  constants,  the 
root-mean-square  roughness  height, respectively ,  for  each  triangular  surface 

element. 


4.  NTSPEC:  Triangular  ground  surface  flag. 
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*Default  Values  for  Surface  Elements  and  Ground  are: 

1.  ISPGRD:  Focusing  ground  computation  flag. 

2.  ERO,  SH20:  Default  values  of  dielectric  constant  and  roughness  height  which 
are  used  in  those  regions  not  specified  by  previously  defined  rectangular  and 
triangular  areas. 

3.  ERG,  SH2G:  Dielectric  constant  and  RMS  roughness  height  for  ground 
reflection.  These  parameters  are  specified  only  once,  since  they  are  assumed  to 
be  the  same  for  the  ground  surrounding  all  buildings  and  aircraft. 

♦Building  Parameters  are: 

1.  ID:  two-character  building  identification. 

2.  (XL,  YL),  (XR,  YR):  X,  Y-coordinates  of  left-hand  and  right-hand  edge  of 
face  of  each  building. 

3.  HB:  Height  of  building,  relative  to  bottom  edge,  for  each  building. 

4.  HBOT :  Height  of  bottom  edge  above  Z-axis  reference  plane. 

5.  SH2B,  ERBR,  ERBI:  The  RMS  roughness  height  and  real  and  imaginary  relative 
dielectric  constants. 

6.  BTILT:  Tilt  angle  of  building.  This  angle  is  positive  if  building  has 
positive  Y  coordinates  and  tilts  away  from  the  centerline  r  if  building  has 
negative  Y  coordinates  and  tilts  towards  the  centerline.  Otherwise  this  angle 
is  negative. 

7.  GRNDBD:  Differential  height  factor  of  ground  on  paths  to  and  from  the 
hangar,  i.e.,  height  of  ground  beside  the  runway  relative  to  the  runway  height. 
This  factor  is  positive  if  ground  is  above  zero  height  level  and  negative 
otherwise. 

The  parameters  ERG  and  SH2G  specified  in  item  3  under  "Default  Values  for 
Surface  Elements  and  Ground"  are  also  used  to  obtain  scattering  from  the 
building  surfaces. 

♦Aircraft  Parameters  are: 

1.  ID:  two-character  aircraft  identification. 

2.  (XT,  YT),  (XC,  YC):  X,  Y-coordinates  of  cockpit  and  tail  fin  ends  of 
fuselage  centerline  of  each  aircraft. 
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3.  NACTYP:  Aircraft  type,  for  each  aircraft,  e.s.: 


1  =  747 

2  =  707-320B 

3  =  727 

4  =  DC  10 

5  =  C-124 

6  =  Convair  880 

7  =  Hastings  aircraft 

8  =  Water  tower 

9  =  Small  diameter  pipe 

11  =  C-5A 

12  =  C-141 

A  subroutine  ACTYPE  is  called,  using  the  appropriate  aircraft  type,  to  load  the 
following  aircraft  parameters,  which  are  already  stored  in  computer  memory,  into 
a  suitable  storage  area: 

1.  Area  of  wings 

2.  Radius  of  fuselage 

3.  Length  of  fuselage 

4.  Radius  of  curvature  of  tail  fin 

5.  Width  of  tail  fin 

6.  Height  of  tail  fin 

7.  Height  of  center  of  fuselage  above  ground 

4.  ALT:  Altitude  of  each  aircraft  defined  as  the  height  of  fuselage  centerline 
above  the  ground  (Z-axis  reference  plane).  If  aircraft  is  parked  on  the  ground, 
ALT  should  be  set  to  zero  so  the  program  can  recognize  that  a  default  value 
should  be  used  in  computations. 

5.  GRNDAC:  Differential  height  factor  of  ground,  i.e.,  height  above  zero  of 
ground  near  side  of  runway. 

The  parameters  ERG  and  SH2G  specified  in  item  3  under  "Default  Values  for 
Surface  Elements  and  Ground"  are  also  used  to  obtain  ground  reflections  for 
scattering  from  the  fuselage  and  tail  fin. 

*PARAMETERS  USED  IN  COMPUTATION  OF  SHADOWING. 

*Building  Parameters  are: 

1.  ID:  Two-character  building  identification. 

2.  (XL,  YL),  (XR,  YR):  X, Y-coordinates  of  left-hand  and  right-hand  edge  of 
each  shadowing  surface. 

3.  HBS:  Height  of  shadowing  surface  relative  to  bottom  edge. 

4.  HBT:  Height  of  bottom  edge  of  surface  relative  to  Z-axis  reference  plane. 
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*Runway  Hump  Shadowing  Parameters  are: 


HUMPF ,  HUMPM,  HUMPB:  X ,  Y  ,Z-coordinates  for  the  location  of  the  hump  along  the 

runway. 

The  runway  hump  is  assumed  to  extend  from  the  lower  to  the  upper  edge  of  the 
runway.  The  peak  of  the  runway  coincides  with  HUMPB  but  does  not  pass  through 
HUMPF  and  HUMPB  unless  they  are  symmetrical  about  HUMPM. 

*Aircraft  Parameters  are: 

1.  ID:  Two-character  aircraft  identification. 

2.  SHP0S1:  X,  Y,  Z-coordinates  of  center  of  fuselage  of  each  shadowing 
aircraft  at  the  starting  frame  number. 

3.  SHP0S2:  X,  Y,  Z-coordinates  of  center  of  fuselage  of  each  shadowing 
aircraft  at  the  ending  frame  number  (assumes  linearpath). 

4.  SHACTYP:  Aircraft  type  (see  item  3  under  "Aircraft  Parameters"  above). 

5.  SHVEL:  Velocity  of  shadowing  aircraft  between  SHP0S1  and  SHP0S2. 

6.  SHANG:  Pitch  angle,  angle  between  fuselage  centerline  and  the  X-axis 
measured  in  the  X-Z  plane. 
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DESCRIPTION  OF  PATH  FOLLOWING  ERROR  (PFE)  AND  CONTROL  MOTION  NOISE 

CCMN)  FILTER  EQUATIONS 

The  math  model  output  data  are  processed  through  standard  filters  in 
order  to  assess  the  effects  of  the  time  reference  scanning  beam 
(TRSEO  errors  on  actual  aircraft  movements  and  on  control  surface 
and  stic k  mot i on . 

PFE  is  defined  as  the  theoretical  worst  case  deviations  from  a 
preselected  course  of  an  aircraft  following  microwave  landing  system 
(MLS)  guidance  commands.  CMN  is  that  portion  of  the  error  which 
affects  control  surface,  wheel,  column  motion,  and  aircraft 
attitude.  Rate  noise  is  a  measure  of  the  rate  error  in  a  bandoass 
region  which  would  affect  aircraft  guidance  accuracy. 

F'FE’s,  CMN,  and  rate  errors  are  determined  by  passing  the  1 1  me 
records  through  standardised  filters,  described  herein.  The  filter 
characteristics  are  based  on  a  wide  range  of  existing  aircraft 
response  properties  and  are  believed  to  be  adequate  for  am 
foreseeable  aircraft  as  well. 

Only  PFE’ s  and  CMN  errors  are  currently  used  in  the  anal  cn  ri  m :« t  • 
model  output. 
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The  actual  aircraft  perturbations  caused  b  y  MLS  err  ore  are  7. 
function  of  the  guidance  loop  bandwidth-  The  aircraft  per  '  .  L h  i 
or  F’FE  will  be  estimated  using  a  path  following  .-liter.  The 
transfer  function  of  this  filter  is  as  follows: 
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After  applying  a  bilinear  transform  to  this  function  we  obtain  the 
following  digital  filter  difference  equation  which  is  implemented 
l n  the  model : 
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CMN  is  noise  which  causes  control  surface  motion,  column  motion,  and 
wheel  motion,  but  does  not  significantly  affect  aircraft  position. 
The  CMf-r  s  will  be  estimated  using  a  control  motion  filter.  The 
transfer  function  of  this  filter  is  as  follows: 

s  Wa 

H  ( s  >  =  - * - 

5  +•  Wi  s  +  W-. 

After  applying  a  bilinear  transform  to  this  function  wo  obtain  the 
following  digital  filter  di f  f erence  equation  which  is  implement-" 
in  the  model: 
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Input,  at  time  n 
Output  at  time  n 
Sampling  period 
Critical  frequency 
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COMPUTATION  OF  MLS  SYSTEM  COORDINATES 


Midway  runway  22L 


General  information  required:  Enter  appropriate  values. 


Runway  Length  (ft)  = 

Runway  Width  (ft)  = 

Threshold  Crossing  Height  (ft)  = 


Lr  :=  5346  =6102-756 


Wr  :=  150 


ARDH  :=  55 


Glide  Slope  Approach  Angle  (deg)  =  MGPA  :=  3.6 


Final  Approach  Fix  (FAF)  (nm)  = 
MSL  reference  value  (ft  MSL)  = 
Threshold  MSL  (ft)  = 

StopEnd  MSL  (ft)  = 


FAF  :=  6.241587 


Zref  :=  605 


Zth  :=  605 


Zse  :=  619 


Initializations: 


deg  :  = 


Ztch  :=  Zth  +  ARDH 


MGPA  :=  MGPA  deg 


Runway  Profile:  Enter  X  and  Z  values  along  centerline. 

nrpp  :=  7 


XR  :  = 
i 
0 

550 

1375 

2950 

3750 

5325 

6102 


("rpp  -  1) 

ZR  :  = 
i 

619 

619 

613 

510 

608 

606 

605 


Terrain  Profile:  Enter  Y  value  of  EL  offset  as  well  as  X  and  Z  values 

along  the  offset. 


Y  :=  325 
el 


nop  :=  15 
j  : =  0  . . ( nop  -  1 ) 


XE  :  = 

_ 


4050 

4150 

4250 

4350 

4450 

4550 

4650 

4750 

4850 

4950 

5050 

5150 

5250 

5350 

5450 

ZE  :  = 


606.8 

607.0 

606.7 

606.3 

606.0 

606.0 

605.8 

605 

606 

605 

604.9 

605.2 

605.4 

605.1 

605.2 

EL  phase  center  is  3  feet  above  ground.  Therefore 


ZEPC  ZE  +  3 

j  j 


ZEPC  , ZE  ,ZR 


X  M 


quation  of  EL  phase  center  constrained  oy  .-mown  values  at  threshold. 

value  is  variable. 

E  :=  lspline(XE, ZEPC) 

Zepc(X)  :=  interp(E/XE,ZEPC,X) 

Zc(X)  :=  Zepc(X)  +  tan(MGPA) ■ (Lr  -  X) 


X  :=  min(XE) 

Xel  :=  root ( ( Ztch  -  Zc(X)),X)  X  :=  Xel 

el 

Z  :=  Zepc(Xel)  -  Zref 
el 

ELEVATION  COORDINATES 

X  =  4614.725752 

el 

Y  =  325 

el 

Z  =  8.992095 

el 

★  ★★★A*******************************************************-.  ****************** 


Zfp(X)  :=  Z  +  tan (MGPA)  [x  -  X 

el  el 


R  :=  lspline(XR, ZR) 

Zrs(X)  :=  interp(R,XR, ZR,X)  Zrs(l)  =  619.001092 


X  :=  X 
d  el 


Zrs(Lr)  =  605.976471 


Y  :=  0 
d 

Z  :=  Zrs  X  -  Zref 

d  L  el  . 


★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★Ac************ 

DATUM  COORDINATES 

X  =  4614.725752 

d 

Y  =0 

d 

Z  =  1.72063 

d 


x  :=  o 

XFP  :=  FAF' 6076.1  +  Lr  YFP  :=  0 

1  1 


ZFP  :=  Z  +  tan(MGPA)' 
1  el 


XFP  -  X 
1  el 


XFP  :=  root( (Zfp(X)  -  (Zrs(X)  -  Zref  +  8)),X) 
2 

YFP  :=  0 
2 


ZFP  :=  Z 


el 


+  tan(MGPA) 


XFP  -  X 


el 


0  . . (nrpp  -  1) 


$|XFP  -  XR  j 
i 

i  :=  npp, (npp  -  1)  . . 1 

XFP  :=  XR 

3+-npp-i  i-1 

i  : =  3  . . npp  +  2 

YFP  :=  0 
i 

i  :=  npp, (npp  -  1)  ..1 

ZFP  :=  ZR  +  8  -  Zref 

3+npp-i  i-1 

******************************************************************************* 
FLIGHTPATH  PROFILE 


D-4 


